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Abstract 

> ' 

■ We consider gravitational self interaction in the lowest approximation and assume 

that graviton interacts with gravitational energy-momentum tensor in the same way as 
it interacts with particles. We show that, using gravitational vertex with a preferred 
gravitational energy-momentum tensor, it is possible to obtain a metric necessary for 
explaining perihelion precession. The preferred gravitational energy-momentum tensor 
gives positive gravitational energy density of Newtonian center. We show also that, 
employing "improvement" technique, any gravitational energy-momentum tensor can 
be made suitable for using in gravitational wave equation for obtaining metric which 
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(3JT), explains perihelion precession. Yet the "improvement" leads to negative gravitational 

energy density of the Newtonian center. 

X. 

1 Introduction 

The field theoretical approach to gravitation gives an additional insight into gravity both 
within general relativity and beyond. For example, this approach can help to choose the 
preferred coordinate system, which enable as to see directly from g^ v how the gravitational 
field affects the Cartesian latticework of rigid rods, existed before switching on gravitational 
field. The field theoretical approach suggests that in a weak gravitational field of a Newtonian 
center the influence on rods is isotropic [Thirring (1961)]. The same conclusion follows 
from the dimensional analysis in [Dehnen, Honl, Westpfahl (I960)]. On the other hand, 
the standard Schwarzschild coordinate system suggests that only radial lengths are affected 
by spherically symmetric gravitational field, cf. for example, [Cohn (1969)], see also §22 in 
[Einstein (1916)]. Yet, in the framework of general relativity we cannot get the answer to this 
question without performing first the complicated procedure of identification coordinates x 
in g lll ,{x) with points in the laboratory, cf. §22 in [Misner, Thorne, Wheeler (1973)]. 

If the influence of gravitational field is indeed anisotropic, then the atom in gravitational 
field should lose its spherical symmetry and this must affect its spectrum, cf. [Thirring 
(1961)]. We note also that the ability to define the preferred system can also help us to 
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resolve the problem of nonuniqueness of gravitational energy-momentum tensor [Nikishov 
(2003)]. 

It is commonly accepted that using field theoretical methods in building up gravity theory 
one almost inevitably end up with general relativity, see, for example, [Wald (1986)]. We 
show that at least in the lowest nonlinear approximation this is not the case. One appealing 
way to get a version of gravity theory is to use the S— matrix formalism with propagators 
and vertices without any recourse to gravitational wave equation. This way we are free from 
the assumption that there are purely gravitational Lagrangian and matter Lagrangian; they 
may be inextricable. 

We consider the 3-graviton vertex build with the help of a preferred gravitational energy- 
momentum tensor which yields positive gravitational energy density of the Newtonian center. 
Using this vertex it is possible to obtain metric necessary for explanation of perihelion pre- 
cession. 

It seems natural to assume that gravitational field should have positive energy density. 
In this case the existence of black holes should be impossible. Some arguments in favor of 
this are given in [Baryshev (1999)]. It is likely that in the future the gravitational theory will 
be able to define uniquely the gravitational energy-momentum tensor. 

In my previous paper [Nikishov (1999)] I have considered several gravitational energy- 
momentum tensors. Only one of them gives positive gravitational energy density of the 
Newtonian center. In this paper I use this tensor in 3-graviton vertex. 

In Section 2 a notation is introduced for the building blocks of gravitational energy- 
momentum tensors. It facilitates the comparison of different tensors and different theories. 
In Section 3 we obtain the necessary g pu from our vertex. In Section 4 we show how to get 
different versions of gravitational wave equations by using " improvement" technique for grav- 
itational energy-momentum tensors. In Section 5 we compare the Feynman's gravitational 
Lagrangian [Feynman Moringo, Wagner (1995)] with that of general relativity. 

2 Notation and some preliminary relations 

We define 

= Vi* + hvw, V = h iw ~ -ffwK h = h a a , r]^ = diag(-l, 1, 1, 1). (1) 

If not said otherwise Latin and Greek indices run from to 3. As we are going to compare 
different gravitational energy-momentum tensors, it is helpful to introduce special notation 
for the building blocks of these tensors. We denote 

t > k = r] jk h a ^h^' a ; t i k = 7] jk h a ^h al3 ~<; t i k = r] jk h p h pa , CT r > k = rf k h p h' p ; 
t i k = h a ^h af3 ' k ; t i k = h ja 'Ph k p , a ; r ^ = h ja ^h k a ^ r > k = ^-{h ja ' k + h ka > j )h a ; 

f i k = h jk ' a h a ; t jk = )^{ti\M + h ka ,M); r jh = h' j h' k ; f jk = h jk ' a h aa ' a ; 
r jk = \{h^h^ k + h^Kp'i)- T> k = h jtT , a h ka , a ; r jh = \{h^ k + h k ^)h^ a - 

Y'^rph^h-y, h ti = ^. (2) 



Similarly for terms with second derivative: 

t 3k = T? k h/h; t 3k = r] jk h a p> al3 h; r 3k = V jk h> al3 h a p; r 3k = r] jk h a P' a a h a p] 

t jk = h jk ^h af3 ; r 3k = ^(h ja > k(3 + h ka ' i/3 )h a p; r jk = h a ^ jk h a p] r jk = h? k h/\ 
r jk = h? k K^- t jk = \{h 3a h k + h' ka hj); r 3k = l -{h^\h a k + h k °>\h a 3 ); 

£ jk = ^h ja ' ff a h a k + h ka '%h a j ); t 3k = X -(h^\h k + F^V). (3) 

1 16 a q 

Tensors T 3 , • • • T 3 , and T 3 , • • • T J are obtained from (1) and (2) by substitutions 
h^v — > h^. There are relations 

j- jk = « jfc. ^ jk = I « * jfc. ^ jk = I « c jfc _ ^ Jfe = d j7c . 

rrjk _ 1 » jfe_I & jfe_I = jfe , « jfe. rfjkfjk. A jk _\ « jfc_ | jfc. £ jk = \ ( jk_ It jk. 

1 4 2 2 ' ' 2 ' 2 

4-ik = \iik_\Zik_\9 jk * jk. rl jk = 1 / jfc _ 1 £ jfc , J jfc _ 1 « jfc. 

4 2 2 ' 4 2 2' 

/f i fc - J i fc - n-i k = -t i k - t i k - ■T i k = - t 3k -- t 3k -- r t 
i ' 7 2 ' 4 2 2 ' 

2j* = \l jk_ « jfc. ° jfc = 1 / jk _ 1 h jk _ 1 n jfc o 

2 4 2 2 ' 

Jf-jh = \^jh_\9j k _ll jk V jk. ± jk = 1 / jfc _ 1 * jfc _ 1 " jfc , 9 jfc / 4 n 

4 2 2 ' 4 2 2 ' v; 

and these 

q~ i k =t jk — t jk + - r T jfc =r T" ■ ?fe = — r jk + — r ]k - T jfc =/ r 
4 2 

jfe = | jfe. <j jk = ^ jk_ 8 jfe + l I. 1 jfc. j jk = jk + 7 9 jfc _ ^ jfc = _ § jfc + 1 11 jfc. 

^ ifc = I ^ 9 ,- fc . . fe = _ io jk + 1 11 jfc . n jfe = n jk _ 

^- 2 ^ = _i|ifc + Iii*_Ii jk + ™jk. £jk = __\*jk_ \ ipj* , I ii , i? j*. 

2 4 2 ' 2 2 4' 

_ _ fjk I 11 jfc , 14 jfc- 15 jfc 10^ 1 11 . fe 15 jfe _ 

i t- 4 r t- r , i 2 2 4 ' 

ri fc = -r^ + -T^T^. (5) 
4 

The reversed relations are obtained from (4) and (5) by substitutions T <-> r. 



Next we introduce the Newtonian potential 0. For one Newtonian center = —^r-, for 
several centers 



a 



For these centers the field theoretical approach gives in linear approximation [Thirring (1961)] 

h$ = -4#^Ao, h$ = -2<j>S^, = = -40 = -P>, Vllu = diag(-l, 1, 1, 1). 

(7) 

The Schwarzschild metric in harmonic and isotropic (but not in standard) frames gives (7) 
in linear approximation. The same expressions is obtained by a heuristic method [Dehnen, 
Honl, and Westpfahl (I960)]. Using (7), for the Newtonian centers we find for T jfc 

1 ., 3 ., 6 ., 8 ., 10 ., 12 ., 13 ., 14 ., 15 ., 16 ., 

q- ]k _rj- ]k _rj- jk _rj- jk _rj- jk _rj- jk _rj- jk _rj- jk _rj- jk _rj- jk _ q. 

T* =T jk = W fe (V0) 2 ; T jk =T jk = 160 J -0, fe ; T jk =T jk = -16(V0) 2 5, o 4 ; o; 
f > k =T 3k = W fc 00,«; f Jk =T ° k =f Jk = -1600, l4 ^o4 ; o; T lk =T lk = 1600, jfe ; 



b ., c ., e ., h ., i .. j ., m ., n ., o ., q ., 
j- jk _j- ]k _rj- jk _j- ]k _rj- jk _rj~ jk jk =/ j- jk _rj- jk _j- jk _ g /g\ 



and for r jfc 



* jfc = I 2 jfc = I | jk = i ^ ^ = 7 jfc = 1_6 ^ = , fc 2 _ 5 ^ = 6 ^ 

4 2 4 v y ' 



2 | jfc = 


2 r ]k = 


^Jk = 




1 a 
- T 

2 


jk = £ jk 


=r i k : 


1 d • 
= - r J 
2 


2 ?ik z 


= 85 ifc 00, 


1 

Mi T\ 


r J fc =• 



4 r J* = 4 r = 160 J 0, fc ; r^ k = 2 r >* = 8<5 jfc (V0) 2 ; 



, = 2 j. jk = 1 * = ™ = 2 ° jfc = 2 | ^ = g^^_ (g) 

j, A; = 0, 1, 2, 3. It is assumed in (8) and (9) that r]j k = diag(— 1, 1, 1, 1). Otherwise one mast 
replace r] jk by -r] jk . 

3 Gravitational energy-momentum tensor 

In my previous paper [Nikishov (1999)] I have considered several gravitational energy- 
momentum tensors in the lowest (e.i. h 2 ) approximation. One of them (corresponding to 
general relativity) is obtainable from the Lagrangian 

Thir 1 1 _ _ . 1 - - - - , 



the other one from 



MTW 1 1- - 1 _ 

L = ^[-2^> + + a] " (10) 



(The Lagrangian (10) is eq. (7.8b) in [Misner, Thorne, Wheeler (1973)]; earlier it was used 
by Feynman, see eqs. (3.6.1) and (3.6.5) in [Feynman, Moringo, Wagner (1995)].) These 
Lagrangians are related as follows 

Thir MTW 1 _ _, a _ 

L = L +^[(V^),a-(V^,a)'1. 

Symmetrizing the canonical energy- momentum tensor, obtained from (10) by Belinfante 
method and adding the interaction tensor 

int ., M . 

T 3k =T Jn h n k , we obtain [Nikishov (1999)] 

fi k — \ j \ j- i k + }_ j- q- i k _ i fj- i k _|_ 2 x jfe — 2 T ~* k — 2 T 
327rG 2 4 2 

2 7- 5 i k - t jk + 2 T jk - 2 T jk + 2 f jk - 2 f jk ] + -(f Jn /i„ fc + f fcn /i„ J ). (11) 
Here T J is energy-momentum tensor of particles 

M j 

T ik = V m a u j u k — S(x - x a (t)), u» = dx^/ds. (12) 

a 

The remarkable feature of tensor (11), called in [Nikishov (1999)] MTW-tensor, is that it 
gives positive energy density of gravitational field of the Newtonian center. Moreover, the 
expression for gravitational energy density has the same form as electromagnetic energy 
density of a Coulomb center: 

M i 

t oo = 2T oo 0+ ( V 0) 2 . (13) 

07TG 

Together with (12) we find for slowly moving particles 

M M i 

T oo +t oo = Y^mJ{x - x a (t))+ T 00 + ^(V0) 2 . (14) 



Here we have used eq.(29) below and eq.(7) according to which h 00 = —20. Now we consider 

M 

only two particles m a and m b and drop the self interaction terms in T 00 

r _^ ir Gm a Gm b 

-\m a b{r - r a ) + m b b{r - r b )\[— — — + 



r-rJ r-r b 



Gm a m m _ Gm a m m 

^-d(r - r a ) - — — o(r-r 6 ). (15) 



|r a -r b | |r a -r b 

Then, the probing particle m a , slowly moving in the gravitational field of m b , has the energy 

see terms proportional to 5{f—r a ) in (14) and (15). This is what is needed 

for explaining the red shift according to §4 in Ch 2 in [Schwinger (1970)]. Indeed, let Woo be 
the frequency of an atom outside the gravitational field. Then Woo(l + 0(r a )) is the frequency 
in the field in terms of (coordinate) time t. In terms of observable time t obs = \/\goo\t at 
point r the observable frequency is w^l + <f>(r a ))(l — <f>(r)). In particular, near the point 



of emission r = r a the observable frequency is tUoo, i.e. the same as without gravitational 
field. This is in agreement with Thirring (1961), who considered electron-proton system in 
gravitational field, using observable length and t obs . 

We see that tensor (11) may be the real gravitational energy-momentum tensor and all 
other tensors are only effective tensors, leading to correct expressions for special cases such 
as total energy of a system, or energy density of gravitational wave. 

In view of this properties we try to use MTW-tensor in three graviton vertex. The 
immediate aim is to obtain goo for the Newtonian center in the form 



(?oo = -(l + 20 + 20 2 ) 



(16) 



necessary for explaining the perihelion precession of Mercury, if the linear approximation for 
9iiv = T)ixv + h^ v is given by (7). So we use the gravitation interaction Lagrangian in the form 



int 

L~- 



-h t^ u 
2 * h \iv l 



(17) 



with t^ u given in (11). 

The corresponding general relativity Lagrangian in terms of t^ v is (see [Iwasaki (1971)]) 



13 

E' 

s=l 



(18) 



Here a s are given in Table 1 



Table 1 



s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


a s 


1 


l 

2 


-1 


i 

2 


1 


-2 


2 


2 


-2 


2 


-1 


2 


-4 



(In [Iwasaki (1971)] an = —2; this must be a mistake.) 
If we rewrite (16) in the form 



3oo = -[l+2(0 + c0 2 )], 



(19) 



then we can find from Iwasaki paper the contribution c s = a s R s to the sum c = c s for 
each s = 1, 2, • • • , 13. Hence we can find R s and use them in different versions of theory. So 
we get R s given in Table 2. 



Table 2 



s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


R s 


l 

2 


-2 


-1 


-2 


-1 


l 

4 


l 

2 


i 

2 


-2 


l 

2 


-1 


-1 


l 

4 


l 

4 


l 

4 


i 

2 


e s 


1 


l 

2 


-1 


i 

2 


1 


2 





1 


1 


3 


-2 


-2 


-2 


-2 









The table is enlarged because in general s = 14, 15, 16 also occur. We note that R s is 

determined by Newtonian value of r ^ . For example, from r ^ =r v v =r v v in (9) it follows 
R 1 = R 7 = R 16 and so on. 

To use Table 2 for the case with three graviton vertex (17), we rewrite P k in (11) in terms 
of r. With the help of (4) and (5) we find 



1 



32ttG" 



1 2 

-- r 
2 



jk 



1 4 jk 
+ - T Jk 

2 



+ T 



jk_ 



2 I? jk _ i? ^ + c j* _ 2 e jk + 2 j jfe _ |, fe + 2 ™^_2r + - (T Jri /i„ fc + T fen V) • (20) 

2 

Next, we convert the terms with second derivatives into terms with first derivatives. For 
the term r ^ v we have 

V r = h^rTh^h = (h^hh)^ - h>°2hh a . (21) 

_ a 4 a 4 

Dropping divergence, we write h^ u r ^ v — > —2h^ v t or symbolically r ^ — > —2 r ^ and 
similarly for other terms. Thus in cubic terms we may make the substitutions 

j iiv ^ _2 ^ A 4 ". £ V u > _2 -f ^ u ■ T V > — T ^ v — 7" 1 • T ^ v > — T ^ v — t ^ v ■ 

j \iv ^ _ 1 fiu _ ^ fj,u. £ jiv ^ _ fiv _ 11 \w . 9. jiv ^ _ ^ fti/_ ^ iiv . h fiu ^ ^ jxv ^ jiv . 

j S 1 2 -j 1^1^ k ^12 / Q 

r P<- v _ -J- r A"'. -f A 4 " _ j A tl/ _ ^ A"^. j- ^ u ^ j- A tl/ _ r A 1 ^. r A 1 ^ _2 7- A 4 ". 

™ a^ _> _2 r ^; r M ^ — > - r r ^; r M ^ — > - r r ^; r ^ ^ -2 r ^; 

« ^ _, _ I a^_ ^ 3 ^ or T ^ — > — 7^ r ^. (22) 

Noting that together with r MJy — > — r ^ v — t ^ we also have r Miy — > — r r we find 

"V *-V (23) 

( We note by the way that substitutions r — > T in (22) and (23) give corresponding relations 
for cubic terms expressed in terms of h). 

With the help of (22) and (23) we find that V~ v in (17) can be substituted as follows 

- ^ 6s ^ ^ + V" mhor+ " VV) ' (24) 
i 

e s are given in Table 2. Now it easy to verify that the contribution to c in (19) from the sum 
in (24) is zero: 

16 

J2 e sRs = 0. (25) 
i 

Using the linearized Einstein equation 

-/V™'% + h ja / 1 + /V"V - h' jn + rf n (h a a - h aX ' aX ) = 16nG T ]n , (26) 
we rewrite the interaction term in (24) in the form 

-(T jn K k + T kn hJ) = -^[r 3k - r jk - t 3k + t 3k - t jk + t jk ]. (27) 
2 IQitG 

M 

For the Newtonian centers this interaction term is equal 20 T 00 . 
Now, using (22) we find 

h jk \{T 3n h n k + T kn h n i) ^ ^-^h 3k [-A ^ k + Ah k + 2r^ k -Ar > k + 



2 t 3k + A t 3k -A r 3 % (28) 
These terms contribute 2 to c in (19). Finally we have to take into account that T 3 consist 

M 

of the observable energy- momentum and interaction part. We need here T ■ For it we have 
T 00 = ^m a ^6(£ - 2 a (t)) ~T 00observ (l + ^oo). (29) 

a 

This is because in the presence of gravitational field 

d s 2 = -g 00 dt 2 (l -v 2 ). (30) 

Here v 2 is physical velocity, see §88 in [Landau and Lifshitz (1973)]. Hence 

dx° 1 1.1,, 

ds V-9oo(l-v 2 ) VI -v 2 2 

M M 

The interaction term \ T 00 h 00 = —0 T in (29) together with (27) makes c = 1 in (19). 
Thus, we can obtain the desired value of g 00 using gravitational energy-momentum tensor 
(11) in three graviton vertex (17). 



4 Gravitational wave equation 

If we want to consider the wave equation approach, we can't just add the gravitational energy- 

M 

momentum tensor (11) to T 3 in the wave equation (26). This would lead to an incorrect 
<7oo- We must first "improve" P k by adding some properly chosen expression, which do not 
affect the conservation laws. The simplest possibility is to use (djdk — rjj k d n d n )h af 3h al3 or 
(djd k — rjj k d n d n )h 2 . The first possibility will lead to an expression proportional to 

5 2 d k 

(j- jk_ j- jk_ q- jfe_|_ q- 3 ky digcarci this possibility because it would change the coefficient 

5 

in front of T 3 in (11), which gives the correct energy-momentum tensor of gravitational 
plane wave, cf. §107 in [Landau and Lifshitz (1973)]. 

Using the second possibility, we add to t jk in (11) the expression 

In terms of r this expression has the same form, see (4) and (5). The constant factor on the 
r.h.s. in (32) is chosen so that g 00 has the desired form (16). Thus, we obtain the "improved" 
gravitational energy-momentum tensor 

t l jk P = tjk + tjk- (33) 

Now, for simplicity we consider a single Newtonian center. Outside the matter we find 

_ GM 2 _GM 2 (5 jk 2 Xj x k \ 



and for in (11) 

GM 2 _ GM 2 ( 5 jk 2x jXk] 

Beginning from equation (34) in all equations up to equation (46) j, k = 1,2, 3. From (32) 
(35) it follows 

3GM 2 imv 7GM 2 {5 jk 2x j x k 



We note that tj k , tj k , and differ only by constant factors and 



Sjk_2xjp\ = (3?) 

Thus, the conservation laws t$ v = are satisfied. 

Comparison of t l ™ p in (36) with corresponding Weinberg tensor ( see equations (7.6.3) 
and (7.6.4) in [Weinberg (1972)]) shows that for the considered special case they coincide: 
t l ™ p = t^ u , if (7) is used for calculating t^ u . The expressions for gj k in Hilbert (hP* = 0), 
harmonic and isotropic frames are 

9T = S jk {l - 20) + cf>\h5 jk - ^), gfr = S jk {l - 20) + <?(6 jk + ^), 

g%° = 6 jk (l-2t + ^ 2 ). (38) 
In all these frames goo are the same as in (16). The differences in gj k are proportional to 

A* + JW = 2(($tt-^)) A, = 5, (39) 

i.e. these (7^ are related by gauge transformations. This is in agreement with the fact that 
in considered approximation Weinberg gravitational energy-momentum tensors in harmonic 
and isotropic frames coincide [Nikishov (2003)]. 

Defining as usual = — ^rj^t, we get from (36) 

pmp _ GM 2 imp _ GM 2 (5 jk _ 7x jXk \ 

to ° ~ 47rr 4 ' tjk ~ An \r* J ' ^ 

Using improvement technique we can make any gravitational energy-momentum tensor 
suitable for using in the wave equation. In this way we obtain another version of theory in 
considered approximation. For example, for Thirring (1961) tensor we have instead of (35) 

Th 7GM 2 _ GM 2 (2x ]Xk 5 jk \ 



8vrr 4 ' jk ~~ 8tt \ r e ' ' r\ 
The "improved" tensor has the form 

C pTh = C + 3t, u , (42) 
where is given in (34). From (41), (42) and (34) we find 

.impTh _ 19GM 2 _ impTh _ 23GM 2 ( 2 Xj x k 5 jk 



From here for barred quantities we get 

zimpTh _ GM 2 pmpTh _ GM 2 ( 23XjXk &jk_ 
f 00 - 47rr 4 > T jk ~ A7r [ r S + r 4 

As seen from (40) and (44) ? Q ™ pTh = t^ v '. Hence, goo are the same in both versions, but gj k 
are different. Using the wave equation in the form 

- K»/ + K* a v + hua% = 167TG7V, (45) 
we find in Hilbert gauge (h tiV \ v = 0) 

<$ = 5 jk (l - 20) + G 2 M 2 + H^*) . (46) 

In conclusion of this Section we note that the nonuniqueness of gravitational energy- 
momentum tensor in general relativity is connected with the nonphysical degrees of freedom 
taking part in its formation [Nikishov (2003)]. If we could define the preferred frame as the 
one formed by only physical degrees of freedom, then the problem of uniqueness could be 
solved. It seems reasonable to assume that for the Newtonian center in linear approximation 
the preferred metric is given by (7) and in the second approximation it may be the harmonic 
one. 




5 Remark on Feynman gravitational Lagrangian 

The Feynman approach to building up gravity theory seems to me quite natural because he 
uses throughout h^ u as gravitational variable, not some function of it. His cubic term of the 
Lagrangian (see equation (6.1.13) in [Feynman, Moringo, Wagner (1995)]) in our notation 
has the form 

o 1 -L.l a 1 i lc Iff i 3 I 1 m n o p 

= ~ 64^G * 8 T jk + 4 T jk ~ 2 T jk ~2 T jk+ T jk ~4 T jk+ 2~ T jk+ T i k ~ 2r i k+ T i k 

,3 1 4 10 1 11 14 16 

+ (r Jk - - r jk - 2 r jk + - r jk + 2 r jh - r jk )}. (47) 

This expression have to be compared with cubic term in the expansion of general relativity 
Lagrangian — 16 ^ G ^/—gR. My calculation gives the following result for this term 

I jk 1 a lb 1 c If 1 g h i j 1 I 1 m n 

~ 16nG 8 T jfc _ 8 T 4 T jfe _ 4 T ~* k ~ 2 T T T T jfc ~~ 4 T ^ 4 T ^ k ~^ T jfe 

1 o p 3 9 ll 32 1 3 1 4 35 1 6 

— 2 r r J fc — 2 T jfc 4 T _ 8 T jfc _ 2 T jfc 8 T jfe 4 T jfc _ 2 T 

3 7 8 1 9 10 1 11 12 13 14 15 1 16 . . 

2 T i k+ T ik - ~T jk+ T jk~ ~T jk+ T jk- T jk- T jk-^Tjk + -T jk)\ (48) 

Using reductions to first derivatives (22) and equation (23), we get for (48) the expression 

13 

(18). For Feynman's Lagrangian (47) the reduction gives (18) with substitution 4 r j k — > 

2 t" 3 jk + 2 r jfc.Thus, it seems that Feynman's Lagrangian (47) differs essentially from (48). 
Evidently Feynman do not think so. In any case his Lagrangian, used as a vertex, explain 

13 14 

perihelion precession. This is seen from the fact that for Newtonian centers r jk =t jk, see 
(9) and text below Table 2. 



6 Conclusion 



Our considerations show that the ways are open in search for gravitation theory with unique 
gravitational energy-momentum tensor, giving the positive energy density. It is still not 
clear which coordinate condition better exclude the nonphysical degrees of freedom; is it 
harmonic, Hilbert or some other condition? Will the future gravitational theory be the theory 
of gravitons (spin-2), or gravitons and spin-0 field (as in [Baryshev (1999)] and [Logunov 
(1998)])? 
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